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, Abstract. Experimental verification of the existence of gravimagnetic fields 

generated by currents of matter is important for a complete understanding and 
|0 ' formulation of gravitational physics. Although the rotational (mirinsic) gravimagnetic 

field has been extensively studied and is now being measured by the Gravity Probe B, 



the extrinsic gravimagnetic field generated by the translational current of matter is less 
well studied. The present paper uses the post-Newtonian parametrized Einstein and 
^j^^ light geodesies equations to show that the extrinsic gravimagnetic field generated by 

the translational current of matter can be measured by observing the relativistic time 
delay and/or light deflection caused by the moving mass. We prove that the extrinsic 
' gravimagnetic field is generated by the relativistic effect of the aberration of the gravity 

■ force caused by the Lorentz transformation of the metric tensor and the Levi-Civita 

connection. We show that the Lorentz transformation of the gravity field variables is 
equivalent to the technique of the retarded Lienard-Wiechert gravitational potentials 
predicting that a light particle is deflected by gravitational fleld of a moving body from 
its retarded position so that both general-relativistic phenomena - the aberration and 
the retardation of gravity - are tightly connected and observing the aberration of 
gravity proves that gravity has a causal nature. We explain in this framework the 
2002 deflection experiment of a quasar by Jupiter where the aberration of gravity from 
its orbital motion was measured with accuracy 20%. We describe a theory of VLBI 
experiment to measure the gravitational deflection of radio waves from a quasar by 
the Sun, as viewed by a moving observer from the geocentric frame, to improve the 
measurement accuracy of the aberration of gravity to a few percent. 



PACS numbers: 04.20.-q, 04.80. Cc 
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1. Introduction 

A gravimagnetic field, according to Einstein's theory of general relativity, arises from 
moving matter (mass current) just as tlie magnetic field arises in Maxwell's theory from 
moving charge (electric current). The weak- field linearized theory of general relativity 
unveils a mathematical structure comparable to the Maxwell equations [HEllS]. Hence, 
this weak-field approximation splits gravitation into components similar to the electric 
and magnetic field. In the case of the gravitational field, the source is the mass of the 
body, whereas in the case of the electromagnetic field, the source is the charge of the 
particle. Moving the charge particle creates a magnetic field according to Ampere's law. 
Analogously, moving the mass creates a mass current which generates a gravimagnetic 
field according to Einstein's general relativity. Ampere-like induction of a gravimagnetic 
field (gravimagnetic induction) in general relativity has been a matter of theoretical 
study since the Lense-Thirring paper [U O [6]. Now, this problem can be tackled 
experiment ally. 

There are two types of mass currents in gravity [7]. The first type is produced by 
the intrinsic rotation of matter around body's center of mass. It generates an intrinsic 
gravimagnetic field tightly associated with body's angular momentum (spin) and most 
research in gravimagnetism has been focused on the discussion of its various properties 

Our recent publications [71 El El [iD] as well as papers by other researchers [TT] [26] 
give a comprehensive review of various aspects of the intrinsic gravimagnetism. It is 
interesting to note that the intrinsic gravimagnetic field can be associated with the 
holonomy invariance group [27]. Some authors [281 129] have proposed to measure the 
intrinsic gravimagnetic field by observing quantum effects of coupling of fermion's spin 
with the angular momentum of the Earth. It might be worthwhile to explore association 
of the intrinsic gravimagnetism with the classic Hannay precession phase [30] 0. 

The first classic experiment to test the intrinsic gravimagnetic effect of the rotating 
Earth has been carried out by observing LAGEOS in combination with other geodetic 
satellites [33l [34] (see also [35l |36l ET]) which verified its existence as predicted by 
Einstein's general relativity. Independent experimental measurement of the intrinsic 
gravimagnetic field of the rotating Earth is currently under way by the Gravity Probe 
B mission [381 EH] that is expected to increment the accuracy of the preceding LAGEOS 
measurement. 

The second type of the mass current is caused by translational motion of matter. 
It generates an extrinsic gravimagnetic field that depends on the frame of reference 
of observer and can be completely eliminated in the rest frame of the matter. This 
property of the extrinsic gravimagnetic field is a direct consequence of the Lorentz 
invariance of Einstein's gravity field equations for an isolated astronomical system 
[iOj m] embedded to the asymptotically-flat space-time and its experimental testing 
is as important as that of the intrinsic gravimagnetic field. The point is that both the 

f This question is intriguing but nothing has been done to clarify this issue except of some papers by 
SpaUicci [311 [32]. 
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intrinsic and the extrinsic gravimagnetic fields obey the same equations and, therefore, 
their measurements essentially complement each other [7]. Furthermore, detection of 
the extrinsic gravimagnetic field probes the Lorentz invariance of the gravitational field 
which determines the gravity null cone (domain of causal influence) on which the gravity 
force propagates. Experimental verification of the Lorentz invariance (causality) of 
gravity is important for the theory of braneworlds [12] and for setting other, more 
stringent limitations on vector-tensor theories of gravity ^3] . 

Ciufolini [H] proposed to distinguish the intrinsic and extrinsic gravimagnetic fields 
by making use of two scalar invariants of the curvature tensor 

h = RaPfivR'^^^'^ , (1) 
1 2 = RaP^^R'^'^'^'' E^jjpa^ , (2) 

where Ra/Sfiu is the curvature tensor, E^i,p„ is the fully anti-symmetric Levi-Civita tensor 
with -E0123 = +\/—g, and g = det{g^^) < is the determinant of the metric tensor. 
Ciufolini [H] notices that weak gravitational field of an isolated astronomical system 
yields /2 = if the intrinsic gravimagnetic field is absent. However, making use of 
I2 is just one of many possibilities to single out the intrinsic gravimagnetic field. This 
question deserves further, more detailed study probably in the spirit of Petrov's algebraic 
classification [15] of gravitational fields or by making use of other techniques (see, for 
example, [T6l 1^ ). One should not confuse the invariant I2 with the gravimagnetic field 
itself. The gravimagnetic field is generated by any current of matter. Hence, I2 = 
does not mean that any gravimagnetic field is absent as assumed by some researchers 
[17] . Equality I2 = only implies that the gravimagnetic field is of the extrinsic origin 
(/i 7^ 0), that is generated by a translational motion of matter. 

The goal of the present paper is to show that the extrinsic gravimagnetic field can 
be measured in high-precision relativistic time-delay experiments conducted in the Solar 
system where light (photon, radio wave) interacts with a moving gravitational field of 
a massive body [18]. The relativistic light deflection in the rest (static) frame of the 
light-ray deflecting body is well-known and was calculated in optics by Einstein |49j . 
Shapiro [50j has derived the relativistic time delay for radio waves, and the general 
relativistic problem of the gravitational deflection and delay of light by an arbitrary 
moving point-like and spinning mass has been solved in our papers [HI El]. Light-ray 
deflection experiments conducted in the gravitational field of a moving body allow us 
to study the global Lorentz transformation properties of the gravitational field since 
the angle of the gravitational light deflection is essentially an integral phenomenon that 
is not reduced to a local experiment in a particular laboratory frame where the effect 
of the gravitational field is reduced to the tidal force. Because both gravity and light 
participate in the gravitational deflection of light, one has to develop a formalism that 
is capable to trace the Lorentz transformation of the gravitational field variables and 
to distinguish it from that for light. The problem is that conventionally, in general 
relativity the speed of gravity is denoted by the same letter as the speed of hght 
making an erroneous impression that gravity is associated with electromagnetism. The 
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fundamental speed c indeed enters many fundamental equations of modern physics but 
it has different physical facets which should be clearly distinguished to avoid confusion 
in the interpretation of physical laws observed experimentally [52]. 

To separate relativistic effects produced by the Lorentz transformation of 
gravitational field from those caused by the Lorentz transformation of electromagnetic 
field, we introduce in section 2 a post-Newtonian parameter e and put it in front of all 
time derivatives of the gravity field variables of Einstein's equations by making use of 
their re-scaling 

d d 

3 ^ ^3 , V ^ ev , > e— , (31 

where j and v are the matter current and velocity respectively. Parameter e can be 
interpreted as a label which goal is to track down the possible deviation of the null 
cone in general relativity from that in electrodynamics that can be characterized as 
the ratio e = c/cg with c being the constant speed of light, and Cg taking a range of 
values from Cg = oo (Newtonian-like gravity) to Cg = c (general relativity) [53l [Slf [55]. 
This parametrization of general relativity preserves the Newtonian limit (e = 0) of 
the Einstein equations and assumes that gravitational potentials and forces obey (in 
the background flat space-time) the Lorentz transformation that is parametrized by the 
speed of gravity Cg. If e were different from unity, the relativistic Lorentz transformation 
for gravity could be different from that in electrodynamics. This difference can be 
tested in the light-ray deflection experiments conducted in the field of a moving 
gravitating body because in such kind of experiments light interacts with time-dependent 
gravitational field and, hence, both speeds of interaction, c and Cg, are involved. The 
plausible difference of the speed of gravity, Cg, from the speed of light, c, might be 
an indicator that the gravity force violates the principle of causality, that is propagate 
information faster than light. We demonstrate how to test this principle in the high- 
precision light-ray deflection experiments in the time-dependent field of a moving body. 

In section 3 we show that the parameter e is indeed a marker of the Lorentz group 
of transformation of the gravitational field variables, given by the matrix A^(e) that 
depends on the speed Cg = c/e. This speed defines the ultimate speed of propagation of 
the gravitational field in Minkowski background space-time. The Lorentz group of the 
electromagnetic field is described by the matrix = Ap{e = 1), and depends on the 
speed of light c in vacuum, so that the effects associated with the null characteristics of 
the electromagnetic and gravitational fields can be clearly discerned by measuring the 
numerical value of the parameter e. 

In section 4 we analyze the propagation of light, using the formalism from the 
previous sections. Utilizing this formalism, we show that relativistic equations associated 
with the propagation of light distinguish between Cg and c already in linear terms of the 
order of v/cg and v/c beyond the static Shapiro time delay. The linear v/cg correction 
is explained by the relativistic effect of the aberration of gravity force associated with 
the gravimagnetic property of the gravitational force propagating with the fundamental 
speed Cg from the moving light-ray deflecting body to the light particle (photon) [56] . 
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Although Einstein assumed that the Lorentz groups for gravitational and 
electromagnetic field are identical [571 EH], this theoretical prediction should be tested 
experimentally. The compatibility of the two groups can be measured by an experiment 
in which both gravitational and electromagnetic fields are observed and transformed 
simultaneously from one inertial frame to another. If the experiment is sufficiently 
accurate, then the difference between the two Lorentz-transformation speeds, c and 
Cg , can be measured. Very Long Baseline Interferometric (VLBI) measurements by 
observing the propagation of radio waves (light) in the gravitational field of a moving 
massive body does have sufficient sensitivity. In section 5 , we discuss the physical 
interpretations of the defiection experiments similar to the VLBI experiment in 2002 
September during which Jupiter passed within 3.7' of a bright quasar and the gravity- 
aberration term (50 /xarcsec) was measured to about 20% accuracy (10 /larcsec) via 
retarded spatial coordinates of Jupiter in the gravitational time delay as observed in 
the barycentric frame [591 EH] • A potentially more accurate experiment to measure the 
aberration (causality) of gravity can be done when the Sun occults a bright radio source 
like 3C279, and in section 6 we discuss theoretical details of the experiment. 



2. Linearized Gravity and Gravimagnetism 



Our Cg-parametrized approach allows us to work out successive approximations to solve 
the parametrized Einstein gravity field equations by iterations with respect to various 
small parameters. We shall use the linearized approximation of the Einstein equations of 
general relativity, expanded with respect to the universal gravitational constant G. The 
linearized approximation is currently sufficient for analysis of gravitational experiments 
within the solar system. Terms of the quadratic order in G may be detected in some 
future missions [HI]. When discussing the astrometric applications, the additional post- 
Newtonian expansion will be done with respect to the parameter characterizing the ratio 
of spatial to temporal derivatives of the metric tensor. For slowly moving bodies and 
test particles this ratio is 



d/dx^ 



C 



while for light-ray particles (for which v 

d/dx'^ 



(4) 



(5) 



d/dx^ 

where parameter e = c/cg. If Cg — > cx) (that is e — > 0), the general relativity collapses 
to the Newtonian-like theory in which gravitational waves have infinite speed so that 
the aberration of gravity force is absent and all other gravimagnetic phenomena also 
vanish because the gravitational interaction is instantaneous. We retain in the post- 
Newtonian expansion of observable astrometric quantities all terms up to the order of 
e^, and neglect terms of the higher order since they are not detectable with the current 
astrometric accuracy. 
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In the linearized gravity the metric tensor perturbation h^js = Qap — ^lap of the 
gravitational field contains only the gravielectric, /iqo = (2/c^)$, and the gravimagnetic, 
hoi = —{4:/c^)Ai, potentials [H El E]. The space-space components of the metric 
tensor perturbation hij = (2/c^)$(5jj, and the higher order terms are neglected. In 
the approximation under consideration the gravielectric, E, and the gravimagnetic, B, 
fields are formally defined as [H [2l [5] 

B = V xA, (7) 

where V = 9/(9j is a spatial gradient and parameter e appears in accordance to equation 

We choose to work in the harmonic gauge [H] imposed on the metric tensor. It 
does not restrict physical applicability and interpretation of our subsequent results for 
the light ray defiection angle and the proper time delay which are directly observable 
and, hence, gauge-invariant quantities. In the linearized approximation the harmonic 
gauge condition is reduced to the Lorentz gauge, imposed on the potentials $ and A, 

-^ + V.A^O. (8) 

In this gauge the linearized Einstein equations can be written down as a system of the 
gravimagnetic field equations [2] 

V-E =4:7rGp, (9) 

V B =0, (10) 

Vx£=-l^, (11) 

c at 

„ „ edE AnGe . , , 

where p and jt are the mass-density and mass-current of the gravitating matter defined 
in terms of the energy-momentum tensor T^,^ of matter as p = T°°/c^, and f = T^^/c 
respectively [571 EH] • 

Equations ([6]), ([8]), ([9]) and f|T2|) are equivalent to wave equations for the gravielectric 
and gravimagnetic potentials 

+ V^^ <l>(t, x)= - AnGpit, x) , (13) 

Cg <9t2 ) Cg 

These equations are of the hyperbolic type and describe the propagation of the 
gravitational field with the speed Cg = c/e. Their physically relevant (causal) solution 
are the retarded potentials taken over volume of the body under consideration 

m,x)=G[ pJfiy}^, (15) 

Ait,x) = ^[ 4^1^, (16) 
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where 



1 , 

s = t \x — y 



(17) 



-g 



is the retarded time due to the finite speed of propagation of gravity from the body to 
the observer (see Fig. [1]). In the case of a point-hke (or spherically-symmetric) massive 
body the point y ^ z, where z is the coordinate of the body's center of mass. In this 
case the retarded time equation f[T7|) is reduced to 



where the body's center-of-mass coordinate z must be calculated at the retarded time 
s. This complicates solution of equation (fT8|) . Only in case of a uniform and rectilinear 
motion of the body, when z{s) is a linear function of time s, equation (fT8|) can be solved 
exactly g53] . 

Solutions f|T5|) . f|T6l) are fully compatible with the matrix A^(e) of the Lorentz 
transformation (aberration) of the gravitational field, which depends on the fundamental 
speed Cg as well. Indeed, equations (fT3!) . (fT^ are invariant with respect to the Lorentz 
transformation of the space-time coordinates, the gravitational potentials and the source 
variables using the matrix A^(e) (see §111 and §IV for more detail). Thus, measuring the 
aberration of gravitational force would allow us to measure its relativistic properties such 
as the strength of gravimagnetic field the speed of gravity propagation [591 [60], stc. 
In the case of the instantaneous speed-of-gravity theory, where Cg = oo, equations (fT3|l . 
IHM are reduced to a single Laplace-type equation for the scalar potential $ and both 
the gravimagnetic potential. A, and the gravimagnetic field, B, vanish. The Lorentz 
invariance of the gravitational field in such theory is totally broken and is not defined. 

The retarded potentials can be expanded in the post-Newtonian series. Formally, 
it is expansion in the Taylor series with respect to the parameter e (that is, with respect 
to the speed of gravity Cg). It yields 



Notice that the post-Newtonian expansion fl20l) of the gravimagnetic potential A 
commences on the term being proportional to the speed of gravity Cg. This is because 
the gravimagnetic field does not exist in the alternative theories of gravity where the 
gravity field propagates instantaneously, like the magnetic field would not exist in a 
hypothetical electromagnetic theory with the speed of light being equal to infinity 0. 

It is important to understand that the effects of the finite speed of gravity 
reveal themselves not only through the retardation effects associated with second time 
derivatives in wave equations ( |T3l) . dill), but also through the first time derivatives 
in equations (l6l)- (lT2l) . It makes sense since the first time derivatives control the 



s = t |£c — z{s)\ , 



(18) 



-g 




(19) 



(20) 



I We keep the speed of light, c, constant 
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gravimagnetic and aberrational effects of tlie gravitational field in the near zone which 
are matched smoothly to the radiative-zone effects of emission and propagation of free 
gravitational waves generated by the isolated system. 

For a point-like massive body the mass current j = pv, and the integrals (fT9l) . (l20l) 
are reduced to 

$ = -^ + 0(e^), A=^$ + 0(e3), (21) 

where r = |r|, r = x — z{t), and z{t) is a coordinate of the body's center of mass in 
the global frame of reference. In section [3] we show that the potential A appears as 
a consequence of the Lorentz transformation of the gravitational field from the static 
frame of the body to the global (barycentric) frame. This establishes a correspondence 
between two ways of calculation of the gravimagnetic potential - directly from the field 
equations in the global frame, and by making use of the Lorentz transformation of the 
potential $ from the static to a moving frame. 



3. Lorentz Invariance and Aberration of Gravity 

The transformation between two reference frames moving with a constant velocity with 
respect to each other, is described by the Lorentz transformation. The gravity field 
equations (fT3l) . (IT^ are invariant with respect to the Lorentz group with the matrix of 
transformation A^(e) having the "standard" form [571 158] 

AVe)=7e-(l-A^)"'/' , (22) 
A°(e) = A^o(e) = -7./?: , (23) 

A^(e)=5^^ + (7e-l)^, (24) 

where the boost parameter /3* = ev^/c = v/cg, and v = (f*) is the relative velocity of 
observer with respect to the gravitating body. The matrix of the inverse transformation 
A^(e) = A^(— e). Moreover, the transformation preserves the form of the Minkowski 
metric 

= A';(e)A-^(6)r^,, . (25) 

The Maxwell equations are invariant with respect to the Lorentz group of 
electrodynamics with the matrix of transformation A^, given by [57t [58] 

X'o = l^{l-f3T'^" , (26) 
A° = A^o = -iP' , (27) 

where the boost parameter /5* = v^/c. The matrix of the inverse transformation 
A(;(/?) = \>i;{-(3). Notice that = A°(e = 1). 

If e 7^ 1, the Lorentz transformation of the gravitational field, given by the matrix 
A^(e), is different from that of the electromagnetic field equations, given by the matrix 



A^- = + (7 - 1) ^ , (28) 
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A^. Physically, this means that gravitational interaction has the speed Cg that differs 
from the speed c for electromagnetic field. In the case of the instantaneous speed-of- 
gravity theory, when e = 0, the transformation (l22l)-( l2^ degenerates to A^(e = 0) = 5"^ 
because the space-time manifold is split in the Newtonian-like absolute time and absolute 
space with a single gravitational potential $ residing on it. 

Lorentz transformation between two inertial frames, and x'", in gravity and 
electromagentic field equations are given by two different equations 

x'" = A2e)x'^ , (29) 



x'" = A^x^ . (30) 



a 

Equation fl29l) transforms coordinates in the gravity field equations f|T3l) . f|T^ while 
equation fl30|) transforms coordinates in the Maxwell equations. It is important to 
observe that in the limiting case of a very slow velocity v the spatial part of the two 
boosts (l29l) and (!30|) is reduced to two Galilean transformations ^57j that are not identical 
because Cg ^ c. Specifically, equation (l29l) yields 

x'' = x'- evH + 0{(3^^) , (31) 

while the slow-velocity limit of equation fl30l) is reduced to 

= x' -vH + OiP^) . (32) 

If one deals only with gravity field equations, the spatial coordinates x* can be rescaled, 
X* ex*, so that the parameter e = c/cg in equation (!3T|) can be eliminated from linear 
terms and the difference between Cg and c will reveal itself only in the post-Newtonian 
terms of order of (f/c)^ and higher 0. Gravitational light-ray deflection experiments 
observe interaction between light and gravity which depends on the spatial distance 
|a; — 2| between coordinate x of the light particle and that z of a massive, light-ray 
deflecting body. Transformation from one to another reference frame changes their 
coordinates differently according to equations (|3Ti) and (1321) . Hence, the parameter e 
can not be eliminated from the Galilean transformation of the relativite distance |a3 — ^| 
by the re-scaling of the spatial coordinates z of the body. Hence, it can be measured 
already in the linear post-Newtonian correction of order of ev/c caused by translational 
motion of the body and appearing explicitly beyond the static part of the Shapiro time 
delay and/or standard Einstein's light-ray deflection equation [621 [63]. 

Since the the current light deflection experiments are not sensitive enough to 
measure the deflection of light caused by the body's spin, we shall neglect the intrinsic 
gravimagnetic fleld associated with the rotation of the body in the remainder of this 
paper. This approximation eliminates from equations (J9l)- (!T^ all rotational currents, 
(see for more discussion of this intrinsic rotation field). The equations fl9|)- f|T^ solved 
in the static frame of the light-ray defiecting body show that in this frame the potentials, 

I Equation (|3ip can be also presented in the Newtonian-like form by replacing time t t — et. If 
e 7^ 1 the dynamic time t — et in gravity is different from atomic time t in electrodynamics. Hence, 
measuring motion of photons in time-dependent gravitational field allows us to find out the presumable 
difference between the two times, r and i, that is to measure e = c/cg (see |55j for more detail). 
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$' = (cV2)/i^o = -GM/r', A' = -{^/A)h'Q^ = 0, and the fields, E' = -V$' ^ 0, 
= V X A' = 0, where r' = |a;'| is the radial coordinate of a field point in the static 
frame, and M is mass of the body. Transformation of the static field potentials to the 
moving frame is obtained by applying the Lorentz transformation to the metric tensor 
(gravitational field) perturbation 

h^p = K{e)k%e)h'^, , (33) 

and is given more explicitly by equations 

$ =^l[{l + ^l)^' + A{(3,-A')\ , (34) 



A = 7,A' + -fl 



$' + (A ■ A') 

7e + 1 



/3e . (35) 



We emphasize that the gravitational potentials are transformed in accordance with the 
law of the Lorentz transformation of the gravitational field fl33l) so that it is the gravity 
null cone which enters these transformations. The light cone enters the law of the 
Lorentz transformation of the electromagnetic field (light). The presence of the gravity 
null-cone terms is traced by the parameter e in all equations, and each physical effect 
directly associated with the finite speed of gravity can be unambiguously distinguished 
from that coupled to the light cone. 

Comparison with the PPN metric [64] shows that our approach formally yields 
the following values of the PPN parameters [5l]: 7 = j3 = 1, ai = 8(e — 1), 
^ = Ci = C2 = Cs = C4 = ^^2 = ^3 = 0. We emphasize, however, that measurement of e is 
not equivalent to an independent measurement of ai in the PPN formalism framework 
[65], since the PPN formalism assumes that e = 1 while ai 7^ 0. This means that 
the two parameters have different physical origin (see next section) and coincide with 
each other accidentally due to the specific of the parameterizations adopted in pi] and 
in the present paper. Besides the measurement of ai is impossible without making 
an assumption about the existence of a global preferred frame, usually associated with 
the isotropy of the cosmic microwave background radiation (CMBR) [65]. Existence 
of such a "preferred frame" in cosmology should not be related to a possible violation 
of the Lorentz invariance of gravitational field. It simply confirms that our universe 
is homogeneous and isotropic on large cosmological scales. In contrast to the PPN 
formalism, the existence of the relativistic effects parametrized by e do not depend 
upon the "preferred-frame" assumption. Evaluation of the numerical value of e itself is 
rendered through the measurement of the retardation of gravity eSect in the gravitational 
time delay of light by moving bodies as shown in [SH [591 EO] and in equations (!671) - (!69ll 
of this paper. 

The Lorentz transformation fl35|) generates the gravimagnetic potential A and the 
gravimagnetic field B associated with the translational motion of the (non-rotating) 
light-ray deflecting body. Neglecting the quadratic terms in equations ( IMl) . ( l35l) one 
obtains in the linear approximation that 

$ = + o (e^) , A = -^v + O (e^) , (36) 

r Cg 
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where r = |r|, r = x — z(t), and z{t) is a coordinate of the body's center of mass in 
the barycentric frame. We emphasize that the Lorentz transformation fl34|) . fl35l) of the 
gravimagnetic potentials assumes that at each instant of time the motion of the body 
is approximated by a straight hne with constant velocity v = dz/dt. We also bring the 
reader's attention that the equations (1361) are fully compatible with, and have the same 
physical content as the post-Newtonian solutions ([T9l) , fl20|) of the wave equations (|T3|) , 

m- 

The Lorentz transformation of the static field E' of the body generates the 
gravimagnetic field 

B = —{vxE) , (37) 

where E = — V$. Equation ( l37|l results from definition and equation ( 136|) . The 
Lorentz transformation (l34l) . (!35|) describes change (aberration) in the arrangement of 
the gravitational field lines of the body, measured in two different frames. Equations 
( !36ll . ( 1371) describe this gravimagnetic effect of the extrinsic gravimagnetic field in the 
linear approximation with respect to the parameter e = c/cg. 

To summarize this section, in classical electrodynamics a uniformly moving charge 
generates magnetic field. This is because the Lorentz transformation generates electric 
current which produces the magnetic field. The resulting magnetic field is real, since it 
couples to electric field in a gauge-invariant way, and can be measured. Its observation 
confirms that electromagnetic field is Lorentz-invariant and its speed of propagation 
is c [66]. Similarly, the gravimagnetic potential ( l36l) and the gravimagnetic field ( 1371) 
are solutions of the Einstein equations and lead to gauge-invariant observable effects 
which can be measured in gravitational experiments. It provides a test of the Lorentz 
invariance of gravitational field and measurement of the fundamental speed of gravity 
Cg = c/e which controls the causal property of gravity (see section 15.61) and the 
gravimagnetic effects. 



4. Lorentz Transformation of Gravity and Light in the Shapiro Time Delay 

4.1. Equation of Light Propagation and Its Lorentz Transformation 

In the general theory of relativity, the equations of motion for light propagating in 
vacuum are null geodesies, and the general relativistic equations of light propagation 
parameterized by coordinate time t are [5^ 

^ = c^fc^r {kT% - r;j , (38) 

where /c^ = (1, k) is a null four-vector {rj^^^k^k" = 0) of a photon, and k = (fc*) is the 
unit Euclidean vector that is tangent to the unperturbed photon trajectory. The null 
vector k'^ is transformed in accordance to the Lorentz group of the transformation of 
electromagnetic field 

= AtA;^ , (39) 
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where and k'^ are components of the null vector of photon in the static and moving 
reference frames, respectively, and A'^ is the matrix inverse with respect to A'^. The 
Levi-Civita connection describes the force of gravitational attraction and is transformed 
in accordance with the Lorentz group of transformation of gravitational field, that is 



If the parameter e 7^ 1, the Lorentz group of electromagnetic field does not coincide with 
that of gravitational field. Hence, the equation fl38l) must be modified in order to keep it 
form-invariant. It is beyond the scope of this paper to derive this generalized equation 
of light propagation in a particular theory of gravity with broken Lorentz invariance of 
gravitational field. This problem has been tackled in our papers [63l [67] . The goal of 
the present paper is to test whether the Lorentz group of transformation of gravitational 
field coincides with that of electromagnetic field. To achieve this goal we assume that 
equation (l38l) is valid in the observer frame with respect to which a light-ray deflecting 
body is moving. Then, we take electromagnetic and gravitational field variables in the 
static frame of the body and apply the Lorentz transformation with e = 1 for light 
and e 7^ 1 for gravity, thus, obtaining their transformed values which are substituted to 
equation fl38l) . The obtained equation is parametrized by parameter e, and measuring 
e in gravitational time-delay experiments allows us to measure the speed of gravity Cg 
with respect to the speed of light. 

Final notice relates to the nature of the light-ray geodesic equation (1381) . Some 
researchers (see [65], [75] and references therein) state that because this equation describes 
propagation of light it can not contain any information about the speed with which 
gravity propagates. Such statement is valid only in case when the light-ray deflecting 
body is at rest both with respect to observer and to the source of light. However, the 
right side of equation ([381) is the gravity force acting on a light particle, and in case when 
the body is moving, the force of gravity does not propagate instantaneously from the 
body to the light particle. This propagation of gravity is in the form of the " acceleration- 
dependent" gravitational field when the light particle moves in a far (radiative) zone 
of the body where transverse-traceless gravitational waves dominate. In the near (non- 
radiative) zone of the body's gravitational field the gravitational interaction propagates 
to the light particle in the form of the "velocity-dependent" gravitational field. The 
acceleration-dependent and the velocity-dependent components of gravitational fields 
always present together in the most general case of accelerated worldline of the light- 
ray deflecting body. In case of a uniformly moving body, the acceleration-dependent 
part of gravitational field of the body is absent and transverse-traceless gravitational 
waves are not emitted. The radiative zone is absent and the light particle always moves 
in the near zone of the body described by its velocity-dependent component. But it 
does not mean that the inetraction of the light particle with the velocity-dependent 
component is instantaneous. The gravitational field still obeys the causal property and 
the light particle can "feel" the gravitational field only after it crosses the gravity null 
cone as shown in Fig. [2] which is based on the result of the integration of the light-ray 




(40) 
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propagation equation (1551) as shown in section 14. 3[ 

Some researchers also states erroneously that the principle of equivalence tells us 
nothing about the speed of gravity [651 [75]. This is because the principle of equivalence 
eliminates gravitational force acting on a test particle at each point of its world line 
by making a local coordinate transformation to a free-falling coordinate system. This 
statement is true, but it is not applicable to the light-ray deflection experiments which 
are essentially non-local and, as such, test the global properties of the gravitational 
field which can not be measured in local experiments. Gravitational field can not be 
eliminated along the entire light-ray trajectory by a single coordinate transformation. 
Therefore, when the body is moving, a light particle propagates through non-zero, time- 
dependent gravitational field which can not interact with the particle instantaneously if 
the speed of gravity is finite. The light particle (photon) propagates through a sequence 
of the gravity null cones depicted in Fig. [2] which describes the interaction of the 
light particle with the gravitational field of the body in the background flat space- 
time. Gravitational light-ray deflection angle and time delay are integrals taken along 
the unperturbed light null cone but the gravity null cone eneters this calculation as 
well through the time derivatives in the Christoffel symbols normalized to the speed of 
gravity Cg (see next section). Our calculations reveal, that numerical values of the time 
delay and the deflection angle of the light particle, after it passes the light-ray deflecting 
body, are determined at each point of the light-ray trajectory (for instance, points 3, 
4, 5 in Fig. [2]) by the retarded positions of the body (positions C, D, E in Fig. [2]) 
connected to the light particle by the gravity null cone. The retarded position of the 
body must not be confused with the point of the closest approach of the light particle 
to the body. The point of the light's closest approach lies on the light null cone and 
can be spatially close to the retarded position of the body at the time of observation. 
However, the equation defining the point of the closest approach does not follow from 
the result of the integration of the light ray in the gravitational field of a moving body, 
and appears nowhere in the mathematical description of the observable quantities like 
the time delay and the total deflection angle of light. 

This mathematical prediction of the retarded position of the body in the expressions 
for the time delay and the deflection angle follows from the causal nature of the gravity 
field equations, and it leads to the inevitable conclusion that at the time of observation, 
t, the photon received by observer (point 5 in Fig. [2]) is deflected and delayed by the 
body's gravitational field from the retarded position (point E in Fig. [2]) of the body, z{s), 
taken on the "last" gravity null cone (see equations fl671) and fl69] below). If the speed 
of gravity Cg is infinite, the body deflects (and delays) light by its gravitational field 
taken at the time t, when it reaches observer, which means that the gravity propagates 
information about position of the body instantaneously, and does not obey the causality 
principle. This consideration reveals that the gravity null cone (determined by the speed 
of gravity) is an essential part of the light-ray deflecting experiments conducted in the 
field of a moving massive body (see section 15.61 and Fig. [3] for further details) which 
space-time structure can be measured in the ultra-precise observation of the deflection of 
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light coming to observer from a star/quasar. This study refutes the arguments presented 
in [651 EH] that are based on misconceptions in understanding of the interaction of 
electromagnetic and gravitational fields. 



4-2. Lorentz Transformation of the Levi-Civita Connection 

The Levi-Civita connection in a static frame is given by the equations 



^lO p/O -p/i 

00 ~ ij 

■^lO -p/t 



rs; = r:^ = r;;, = , (4i) 

(9$' 

^oi - 00 - Q ,i ' \^^) 



+ + ^43) 

where $' = (c^/2)/;,qq = —GM/c'^r' , and r' = \x'\. As follows from equations (j4T|) - 
(|43|) . the Levi-Civita connection is associated in the static frame with the force of 
gravitational attraction [STJEH]: it is not associated with electromagnetism. Therefore, 
its transformation in equation (HO!) from the static to a uniformly moving frame must be 
done with the matrix A^(e) pertained to the Lorentz transformation of the gravitational 
field. Substituting the matrix of the Lorentz transformation fl22l) - fl2^ to equation PUI) 
and making use of equations fHTl)-(H3l) along with equations (IMl) and ( 135|) . yield the 
components of the Levi-Civita connection in the moving frame 

- - (44) 

^ °° " c,dt ' 

where the parameter e = c/cg explicitly appears in front of the time derivatives in 
equations fH4l) . fj46|) - fj48|) as a direct consequence of the Lorentz transformation law fHOj) . 
Its appearance complies with the phenomenological parametrization rule in equation 
([3]), tracking the presence of the fundamental speed of gravity in various gravitational 
equations. 

It is important to compare equations fj44l) - fj49|) with the definition of the Levi-Civita 
connection adopted in the PPN formalism [64J which postulates that the transformation 
of the connection obeys the "electromagnetic" rule 

r^. = K^'.^T^u , (50) 
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where is the matrix of the Lorentz transformation of electromagnetic field. 
Straightforward calculations reveal that 

ro - - 1^ (51) 
r° - - — (52) 

'dA^ dA'\ 19$ 



(9$ 9$ (9$ 

which should be compared with equations flH|) - (H^ . One can see that the PPN-defined 
Levi-Civita connection fl5T]) - fl56l) assumes explicitly that the parameter e = 1, which 
excludes it from the set of the parameters that probes the relativistic transformation 
properties of the gravitational force (the Levi-Civita connection). 

Obviously, our e-parametrization of the Einstein equations [M] will make all time 
derivatives of the gravitational potentials vanish if the speed of gravity Cg ^ oo. This 
is intuitively pleasing since the limit Cg —>■ oo corresponds to the case of the infinite 
speed-of-gravity theory where the gravitational interaction is instantaneous, and both 
the gravity field equations and equations of motion of test particles, which depend on 
the Levi-Civita connection, should not contain any time derivatives of the gravitational 
field. 

In contrast, the PPN formalism [64] demands that the equations of motion of test 
particles depend on the connection (l5T |) -(l56 l) that contains non-vanishing first time 
derivatives of the gravitational field even in the limit Cg — >■ oo because equations (15T|) - 
( I56ll do not depend on the fundamental speed of gravity Cg explicitly. In other words, the 
PPN formalism orders that the gravitational force propagates with the speed of light and 
in this framework, the speed of gravity cannot be tested. We, thus, consider the PPN 
equations (!5T]) - (!56|) essentially incomplete because the degree of freedom corresponding 
to the propagation of gravitational force that may differ from the speed of light is 
missing. 



4-3. Gravitational Time Delay and Deflection of Light 

The double integration of equation ( l38i) along the unperturbed light ray yields the time 
of propagation of light from the point Xg to the point, = x*(t), on the trajectory of 
the light ray 

t-to = -\x-xo\+A{e), (57) 
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where the relativistic time delay 



A(e) 







dr da 




J to J ~oo 



is a function of the coordinates of emission, Xq, and the point, 
particle is located at the time t. Integration in equation 
unperturbed propagation of the photon 

XN{t) = a^o + ck{t - to) , 



(58) 

X*, where the light 
is along a straight line of 



(59) 



where to is time of emission, Xo is the coordinate of the source of light at time to- By 
substituting equations (H^ - P9l) into the time delay equation (l58l) and making use of the 
gauge condition in equation ([H]) in order to replace the time derivative of the potential 
$ to the divergence of the potential A, we obtain the following form 



A(e) = 2 / <l>(r,a;7v(r))rfr-2 



to 



1 

1 - - 

e 



dr 



to 



V-A(a,x) 



da+0 {5^) , (60) 



where 6 = e — 1 and all terms of order 6"^ or higher are omitted. The potentials $ and A 
are given by equation fl36|) and are taken on the unperturbed light-ray trajectory fl59|) : 
that is, X = Xi\f{t). We emphasize that the standard PPN formalism, with e = 1, ignores 
the possible presence of the second term in the right side of equation fl60|) which comes 
from the explicit dependence of FJ^^ on Cg in the right side of equations (l44l) -( l49l) . For 
this reason the PPN formalism is not sensitive to violation of the Lorentz-invariance of 
the gravitational force associated with the transformation of the Levi-Civita connection 
[54] , and is helpless in physical interpretation of the gravitational light-ray deflection 
experiments. 

When light propagates through a gravitational field, it is not only delayed but 
deflected as well. Let us denote a* the angular coordinates of the deflection vector 
referred to the plane of the sky that is orthogonal to the unit vector fc* that defines 
propagation of the unperturbed light ray. In general relativity, the deflection vector a' 
is defined as a spatial derivative of the relativistic perturbation of the electromagnetic 
phase EH] that is proportional to the time delay (1601) . plus the observer- dependent 
relativistic correction taking into account non-integral (local) interference of the light- 
ray particle with the gravitational field perturbations We use this definition in our 
parametrized approach, and define the parametrized value of the deflection angle as 

'9A(e) 



(61) 



where P*-' 
vector /c*. 



6^^ — is operator of projection on the plane of the sky orthogonal to 



§ Exact Lorentz-invariant equation for the total light-ray deflection angle a* in the framework of 
general relativity for the case e = 1 has been derived in our paper j51| and is given by equation (67) in 
there. 
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4.4- Gharacteristic Time Interval and Spatial Domain of the Gravitational Deflection 
of Light 

In the evaluation of equation ( l60l) . the duration of the time integration interval is 
important for better understanding of the gravitational physics of the light-ray deflection 
experiments. The region of space and time interval during which the interaction of light 
particle with gravitating body is relevant can be estimated as follows: A light particle 
is deflected (scattered) by a long-range gravitational force whose potential $ decreases 
basically as 1/r, where r is the distance of the light particle from the light-ray deflecting 
body. The approximate expression for the deflection angle a of a light particle emitted at 
infinity is obtained from equation (iGlll . if one neglects motion of the light ray-deflecting 
body. For an observer located behind the body with respect to the source of light, the 
result is [571158] 

where d is the impact parameter of the light ray with respect to the body that is fixed 
in space at some instant of time ^ Let us assume that the astrometric accuracy for the 
angular measurements is Aa. Then, the region of radius from the body where the 
gravitational interaction of light with the body is significant is determined by solving 
equation 

ttoo - Oi{rs) = Aa , (63) 

where ttoo = 4:GM/c'^d. We assume that d < rg so that we can solve equation fl63l) by 
expanding it in a Taylor series with respect to the parameter d/rg, which gives 

rs^d(^y^\ (64) 



,2Aa 

The time that it takes light to cross this region is 
2r. df2a„Y" 

which gives us a characteristic time of the gravitational interaction of light with the 
body under consideration, assuming an angular tolerance of Aa. 

Current VLBI technology can measure changes in the relative angles between 
sources in the sky with an accuracy ~10 yuarcsec [60], and in the next decade the 
accuracy may approach ~ 1 /iarcsec = 5 x 10^^^ rad [TOl [TH ESj- As an example, 
consider the Sun as a light-ray deflecting body with the light-ray impact parameter 
d = Rq = 7 X 10^° cm - the radius of the Sun. Using a tolerance of Aa = 10 /iarcsec, 
we obtain from equation that Vg — SOOi?© = 2.1 x 10^'^ cm. This means that a 
moving particle of light continuously interacts with the gravitational field of the Sun 
in a region around the Sun with the size which is slightly above one astronomical unit 

II More precisely, the impact parameter of the hght ray is the spatial distance of the closest approach 
of the ray to the light-ray deflecting body. In general, it depends on the coordinate frame used for 
calculations. 
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(1 AU = 1.5 X 10^'^ cm). The light travel time across this region is tg ~ 23 minutes 
which gives the characteristic time of interaction of light particle with the gravitational 
field of the Sun for an assumed astrometric accuracy. Similar calculation for the light 
particle moving in the gravitational field of Jupiter and grazing its limb, shows that the 
characteristic time of light scattering is about tg — 14 seconds, and the region of the 
gravitational interaction of light with Jupiter is about ~ 30 jovian radii. 

These estimates elucidate that the characteristic time of interaction of 
electromagnetic signal and gravity in the light-ray deflection experiments is sufficiently 
large so that this interaction can not be considered as an instantaneous scattering 
by a point-like mass to produce a total deflection angle Ooo- This "instantaneous- 
scattering" point of view about the gravitational deflection of light is commonly used 
in the gravitational lensing theory to map images of the sources of light from the source 
plane to observer [71]. However, it may be a misleading concept for analysis of the 
light-ray deflection experiments in the solar system. 

5. Summary of the Physical Interpretations of the Deflection Experiments 
in the Field of a Moving Gravitational Body 

5.1. Significance of the Velocity- Dependent Corrections Caused by Motion of the Body 

Relativistic corrections caused by the motion of the light-ray deflecting body are directly 
associated with transformation properties of the gravitational field regulated by the 
fundamental speed of gravity Cg which defines the null cone of gravity propagation. For 
example, consider the relativistic defiection of light caused by the Sun and viewed from 
the geocentric frame. The Sun moves in this frame with velocity Vq ~ 30 km/s and 
acceleration Qq ~ 0.6 cm/s^. Let us also assume that position of the Sun z{t) is fixed 
at some instant of time tA, z{tA), and that the total angle of the defiection of light is 
given by equation for a static gravitating body: a = AGMQ/c'^d, where d is the impact 
parameter of the light ray with respect to the body referred to the time tA- 

A light ray grazing the solar limb {d ~ Rq) interacts with the Sun during time 
tg < 23 minutes for an assumed angular accuracy of 10 /iarcsec. Hence, if one ignores 
the geocentric velocity of the Sun it will produce an error 6d = v^tg in the value 
of the impact parameter, which leads to error in calculation of the angle of the 
gravitational light defiection of 5a^ ~ aQiv^tg/ Rq) ~ 105 mas, where a© = 1.75". The 
geocentric acceleration of the Sun, a© = Wq, introduces an error 5d = aQtl/2 to the 
impact parameter and an error 5aa ~ aQ{aQt1/2RQ) ~ 15 /iarcsec to the defiection 
angle 0. These estimates clearly show that, if precision of angular measurements is 
approaching 10 /iarcsec, the relativistic defiection of light and the Shapiro time delay 
in the geocentric frame are affected by the geocentric velocity and acceleration of the 
Sun and, hence, must be taken into account for adequate physical interpretation of the 



^ Proper physical treatment of the acceleration-dependent relativistic effects in the propagation of light 
ray requires more extended mathematical calculations (see, for example, |51j). 
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results of the solar gravity-deflection experiments El- 

Because both light and the gravitating body are moving as the light propagates, 
the instant of time tA associated with the impact parameter d (when the deflection is 
greatest) requires precise evaluation of appropriate time delay and deflection integrals. 
Some claims state [76] that this instant of time is the time of the closest approach 
of photon to the light-ray deflecting body. However, these claims were based rather 
on an intuitive guess than on precise mathematical calculation of t^- The appropriate 
calculation shown below in equations (I66l)-(l69i) for the time delay in the fleld of a moving 
gravitating body (see [SU [591 [63]) shows that it deflects light from the retarded position 
z{tA) = z{s) on its world line where the retarded time s is deflned by equation ([681) for 
null characteristics of the gravitational fleld produced by a moving body. This result is 
in a full agreement with general relativistic prediction that the causality principle for 
gravity is valid and that the gravity propagates and interacts with other particles on 
the null cone. 



5.2. Time Delay and Deflection of Light by a Moving Body 

Integration of equation f[60|) for the time delay can be performed with sufficient accuracy 
under the assumption of uniform and rectilinear motion of the light-ray deflecting body 
with constant velocity v, that is 

z{t) = z{tA) + v{t - tA) + 0{t - U)' , (66) 

where z{tA) is the position of the body at time ^a, with quadratic and acceleration- 
dependent terms neglected. Light propagates through time-dependent gravitational 
fleld as shown in Fig. [2l and the time delay is given by an integral along the 
light-ray trajectory. The evaluation of the integral in the time delay is tedious but 
straightforward, and exhaustive description of this calculation can be found in sections 
3 and 4 of our paper [5l], so they are not repeated here. The result is given by [SH, [55] 

A(e)= -{l+-f)(^l-^k-v^^\n{r~k-r) + 0{6^) , (67) 

where 6 = e—1, r = x — z{s), v = dz{s)/ds is velocity of the body, and the coordinates 
and velocity of the light-ray deflecting body in equation ([67]) are taken at the retarded 
time tA = s where 

s = t - —\x - z(s)\ . (68) 

Equations ([SZD, dSHD were used to draw Fig. [21 It shows that the relativistic time delay 
of photon at each consequtive point on its trajectory is determined by the position x 
of the photon at time t, and the retarded position of the body z{s) on the gravity null 
cone with vertex of the gravity cone on the world line of the body taken at the retarded 

The magnitude of the acceleration-dependent terms in the hght-ray deflection experiments shown in 
[75] is at the thousandths of a picosecond level (< 10~^ yitarcsec). We emphasize that this estimate 
relates to terms of higher order than those taken into account in the present paper. 
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time s in accordance with equation (I68p . Notice that neither the point nor the time of 
the closest approach of the photon to the body appear in the expression fl671) for the 
time delay. 

From now on, we shall treat the point x on the light-ray trajectory as taken at the 
point of observation, and time t as the time of observation (point 5 in Fig. [2]). This 
is because we can not make any measurement of photon before it arrives to observer. 
It is remarkable that the retardation time s of gravity is determined only by the "last" 
gravity null cone connecting the body and the observer at the time of observation of 
the photon t. For this reason, we do not show in Figs. ([3]) and ([3]) the entire sequence 
of the gravity null cones interacting with light like we did in Fig. O The appearance of 
the speed of gravity Cg in equation fl68l) reflects the fact that the gravity force generated 
by the moving body does not interact with light instantaneously but with a finite speed 
because of its causal nature. Had the gravity propagated instantaneously with Cg = cx), 
the time delay equation (167|) would depend on the position of the light-ray deflecting 
body taken at the time of observation t. Graphical presentation of the gravitational 
physics lying behind the time delay equation (|67|) is shown in Fig. [2] and [3] (see also 
graphics in [63]). 

The light-ray deflection vector in the plane of the sky is calculated from equation 
(16T|) and is given by 

r V ■ r 

Cg 

where ^* = P^ r-' is the impact distance of the light ray to the moving body taken at the 
retarded instant of time, and = P^jV^ is the transversal velocity of the body lying in 
the plane of the sky orthogonal to vector k^. For generality, we have added the PPN 
parameter 7 in equations (1671) and (1681) to include possible violations of general relativity, 
caused by presence of hypothetical scalar fields remained after epoch of the cosmological 
inflation [77]. The addition of 7 does not interfere with the Lorentz transformations 
]5T\ ES] or contradicts any of the previous formulations. 

We emphasize that the velocity-dependent effects (i.e. v/cg) appear both explicitly 
in the time delay and light-ray deflection equations (j671) . (j69l) and implicitly through 
the retarded time equation (1^ that displaces the position of the body on its world-hne 
backward in time with respect to the instant of observation t. This occurs because 
equation (jHED is a null characteristic of the gravitational field equations ffT^ . ffT^ 
describing the relativistic effect of the propagation of gravity [59j. The appearance 
of Cg in the time delay equation 0671) is a natural consequence of the fact that the 
Levi-Civita connection fH4|) - fH9|) contains time derivatives which are normalized to the 
speed of gravity Cg = c/e in accordance with the Lorentz-invariant properties of the 
gravitational field that controls the propagation of gravity force both in the near and 
radiative zones in the form of free gravitational waves [631 EH]- In other words, the 
Lorentz invariance of gravity and its finite speed of propagation are physically tightly 
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connected and should not be conceptually separated from each other. 

The physical interpretations associated with the light deflections experiments in 
the gravitational field of moving bodies are also discussed elsewhere (see, for instance, 
[5^ [55| [60l [63]). However, in the remainder of this section we summarize them for 
completeness and for the reader's convenience. 

5. 3. Relativistic Delay Measurement of the Gravimagnetic Field 

Equation fl67|) describes the relativistic time delay of light (radio waves) caused by a 
massive body moving with velocity v with respect to the rest frame of observer located 
at the point x. Gravity is a long-ranged field and the time delay given by equation (l60l) 
is effectively an integral of the gravitational force exerted on the photon along its entire 
trajectory. However, the strongest impact of the time-dependent gravitational field on 
the photon emitted at time to is when the massive body is located in its retarded position 
z{s) taken at the retarded time s, determined by equation ( l68i) . If general relativity 
is correct, then, Cg = c and only the gravielectric potential $ of the body is essential 
for calculation of the time delay (1671) . However, if Cg 7^ c (hence, general relativity is 
not correct), then the gravimagnetic potential A adds a contribution to the delay, as 
shown by equation (l60|) . The second term (the double integral) in the right side of 
this equation is produced by the gravimagnetic field due to the orbital motion of the 
light-ray deflecting body. If general relativity is correct, the contribution of the double 
integral in equation (|60ll to the time delay is identically zero, and the magnitude of the 
gravimagnetic field is given by equations fl36|l and fl37|) with e = 1. Hence, measuring 
parameter e is the test of the presence of the extrinsic gravimagnetic field. This point 
has been misinterpreted by Pascual-Sanchez [IT] who confused the relativistic effect of 
the gravimagnetic field caused by the mass current due to the Lorentz transformation of 
the gravity field variables with the measurement of the classic (non-relativistic) Romer 
delay of light emitted by Jupiter. 

5.4- The Lorentz Invariance of Gravity 

The introduction of parameter e, associated with the Lorentz transformation of 
gravitational field variables, distinguished between the gravitational effects due to the 
fundamental speed of gravity Cg from the special relativistic effects caused by the speed 
of light c. The matrix of the Lorentz transformation of the gravitational field variables, 
A'^(e), depends on the speed of gravity Cg = c/e which is the same constant as that 
entering the parametrized Einstein field equations. If e 7^ 1, the Lorentz invariance of 
gravity is broken with respect to the Lorentz invariance of electromagnetic field and 
contrariwise. Thus, measuring the ratio e = c/cg allows us to test the fundamental 
compatibility of the gravitational field with the causality principle which is postulated 
in general relativity but can be violated in alternative theories of gravity [EU [69] . 
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5.5. The Aberration and the Speed of Gravity 

The e-parameterization of the Einstein equations with the single parameter e helps to 
keep track of any gravitational effect associated with the fundamental speed of gravity 
Cg = c/ e. If Cg = oo, all terms in the gravity field equations depending on time derivatives 
of the gravitational field would be completely supressed and, hence, propagation of the 
gravity force would not be detectable, and the gravimagnetic effect would be completely 
unobservable as well. This interpretation is in a full agreement with the causal (retarded) 
nature of gravity which is revealed in equation (!68|) describing the null characteristic of 
the gravity field equations (fT3|) and ffl^ connecting the field point x with the retarded 
position z{s) of the light-ray deflecting body. The null characteristics of gravity must 
not be physically confused with the null characteristics of light in the light deflection 
experiments. Indeed, light propagates from a source of light to observer while the null 
vector n = r/r, where r = x — z{s), points from the observer, located at x, to the 
gravitating body located at the point z{s) along the null cone of gravity. In general 
relativity Cg = c and the characteristic hypersurfaces of the null cones for gravity and 
light coincide [79]. Nevertheless, null lines of gravity and light propagation can be 
distinguished since they are directed to the source of gravity field (the body) and to the 
source of light which are clearly separated in observations [60l [631 • 

One notices that in some papers the field is considered as a medium in a flat space- 
time with a refractive index different from unity [171 |80l Ell [82]. This interpretation of 
the gravity field should be accepted with a great care since the refraction index of any 
material medium is a scalar and its Lorentz transformation properties are different from 
those for the gravitational field which is a tensor of second rank. Ignoring this subtlety 
can lead to misinterpretation of gravitational effects. 

In light-ray deflection experiments the retarded position of the light-ray deflecting 
body is measured not by observing body itself, but from precise measurement of the 
direction of the gravitational force exerted by the moving body on propagating photons 
in the plane of the sky [5ll [591 [HO]- The retardation of gravity effect can be obtained 
theoretically in two independent ways: 

(i) by making use of the retarded Lienard-Wiechert potentials of equations f[T3l) - f[T4l) 
[5I1E91163]; 

(ii) by solving equations of light propagation in the static frame of the body and, then, 
using the Lorentz transformation, as it is done in this paper and in [83] . 

The two methods lead to the same results given by equations ( [67]) and ( l68i) . The Lorentz 
transformation of the time delay from a static to a moving frame means that both light 
and gravity field variables entering equation (l58ll must be transformed simultaneously. 
The aberration of light transforms vector fc" = (1,A;) of the light ray alone. If the 
aberration of light is taken into account in equation f[58l) by transforming the Levi- 
Civita connection (gravitational force) with Cg ^ c, the time delay can not maintain 
its invariance, and the terms proportional to S = c/cg — 1 will emerge in equations. 
If the Lorentz transformation matrix of the Levi-Civita connection is the same as for 
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light, then, 6 = and the time delay is Lorentz-invariant with the position of the 
light-ray deflecting body taken at the retarded time s given by equation (|68|) with 
Cg = c. Thus, the aberration of light is used in the deflection experiments exclusively 
as a calibrating standard for which speed c is flxed and not measured. The speed of 
gravity Cg is then measured with respect to this standard [62]. It is worth emphasizing 
that all other gravitational experiments rely exclusively upon the JPL [85] or equivalent 
EPM2004 [86] ephemerides of the solar system constructed in the metric system of 
units with a numerical value of the speed of hght c = 299792458 m-s~^ precisely (see 

also [HUEHlEniEolEI]). 

5. 6. The Speed of Gravity, Causality and the Equivalence Principle 

Review article [65] (see also [75] and references therein) argues that the speed of 
propagation of the gravity force from a moving body to the traversing photon is irrelevant 
in the light-ray deflection experiments done in the Solar system, and that the time delay 
(!67|) and the light-ray deflection angle (169|) could be calculated correctly even if the speed 
of gravity Cg would be inflnite. In other words, the author of [65] maintains the point 
that the causal structure of the null cone of the gravity fleld (see Figs. [2] and [3]) is not 
essential to the calculation of the light-ray deflection angle so far as the acceleration 
of the light-ray deflecting body is ignored. This point of view stems from the belief, 
explicitly formulated and adopted in the PPN formalism [61], according to which the 
principle of equivalence demands nothing about the speed of gravity that determines 
the hypersurface of causal influence of gravity. We are told that this is because the 
principle of equivalence operates with derivatives of the metric tensor of the flrst order 
(the Christoffel symbols) while the speed of gravity is determined only by the second 
time derivatives of the metric that appear in the fleld equations of each metric theory 
of gravity (for more detail see section 10.1 in [HI]). Therefore, since the Christoffel 
symbols can be transformed away at any point on the space-time manifold, only the 
second derivatives (tidal forces) can be locally measured, and only these derivatives 
show up speed properties of gravity. 

First of all, we notice that the speed of gravity and the (non-local) causal structure 
of space-time are tightly connected. Gravitational interaction can not propagate 
information about position of a moving massive body faster than the speed of light 
if general relativity is correct. Hence, experimental study of the causal properties of the 
curved space-time manifold answers the question if the speed of gravity has the same 
numerical value as the speed of light or not. The PPN postulate, mentioned above, 
assumes that the causality of gravitational fleld should be relegated only to the gravity 
field equations. However, the causal structure of gravitational fleld is derived not so 
much from the fleld equations themselves but from the exploration of the behavior of 
the set of future directed timelike and/or null geodesies in a given space-time manifold 

l92] . This behavior can be determined, at least in a close neigborhood of any event, 
by the geodesic equations without imposing Einstein's fleld equations [93]. The light 
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geodesies define the causal past of observer, tliat is the region bounded by the past light 
cone in Fig. [31 The causal future of the gravitational field of a massive body, that is the 
region inside the future gravity cone in Fig. [3l is defined by the null geodesies associated 
with the propagation of gravity field from the body. Causal gravitational interaction of 
the body with the observed photon implies that the causal future of the gravitational 
field is not degenerated to a space-like hypersurface that is the photon is gravitationally 
deflected by the body not at the time of its observation but in the past, so that the 
observer and the light-ray deflecting body are connected by a null line which must be 
a solution of the field equations as demonstrated in Figs. [2] and [31 In general relativity, 
where Cg = c, the past light cone of observer must touch the future gravity cone along 
a null line which is a bi-characteristic for the two null cones. 

Second, the equivalence principle tells us that in a local reference frame light 
moves along a straight line [571 EHj- But gravitational light-ray deflection experiment 
is not experiment in a single local frame of reference as photon propagates through a 
continuous sequence of local frames, thus, accumulating the pointwise influence of the 
gravitational field at different parts of the light-ray trajectory. Photon's propagator is 
an integral of the affine connection that is the gravity force shown in the right-hand side 
of Eq. (l38l) . This propagator contracted with the projection operator P^^ , yields the 
integrated deflection angle (169|) which can not be transformed away if the space-time 
is curved, irrespectively of the property of the Christoffel symbols to turn to zero at 
the origin of each local frame. Gravitational deflection of light is a global phenomenon 
which is not reduced to the test of the principle of equaivalence. Notice that even in 
the static gravitational field the principle of equivalence is able to explain only a half 
of the angle of the gravitational deflection of light while another half of the effect is 
due to the curvature of space [571 IMl IM]. In the case when the light-ray deflecting 
body moves with respect to observer the deflection effect gets more sophisticated as 
the magnitude and direction of the deflection angle become functions not only of mass 
but also of the time-dependent position, velocity, etc., of the light-ray deflecting body. 
Therefore, by precise measurement of the angle of the light deflection (or, equivalently, 
the time delay) one can determine how strong the gravitational field of the body is, and 
how fast it interacts with the light ray as it traverses toward observer. If one is able to 
derive position of the massive body from the precise measurement of the deflection angle 
(!69|) and to confirm that the body and the observer are connected by a null line which 
is a characteristic of the future gravity null cone, it gives a direct proof of the causal 
nature of the gravitational field and allows us to measure its finite speed of propagation 
[59] . However, this measurement of the causal property of gravity is impossible if the 
gravitational field is static. Indeed, in case of a static planet the gravitational field 
does not change as photon moves from star to observer, and the causal character of the 
gravitational field is invisible because the body is always at a fixed distance from the 
observer (see Fig. [HA.). 

The situation changes dramatically if the body moves with respect to observer 
because it makes the distance between the body and the observer depending on time 
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(see Fig. IHB). In this case, photon traverses through the gravitational field that changes 
on the light-ray trajectory due to the motion of the body, even if this motion is uniform. 
Had the speed of the propagation of the body's gravitational field been different from 
the speed of light it would have unavoidably led to violation of the causal nature of 
gravity that would be inconsistent with the principle of equivalence as shown in [69] . For 
example, the instantaneous propagation of gravity would imply that one could determine 
current (as opposed to retarded) position of the body on its orbit from observation of the 
gravitational deflection of light, that is the gravitational field would transmit information 
about the body's spatial location to observer faster than the light arriving from the star. 
This violates the principle of causality, and we conclude that correct description of the 
gravitational physics of the light-ray deflection experiment with a moving body requires 
taking into account both the light and gravity null cones as demonstrated in Figs. [1]- 
m and supported by calculations in this paper as well as by the discussion given in 
[54] . Additional arguments backing up the concept of the retardation of gravity in the 
light-ray deflection experiments by a moving body are discussed in our paper [63] both 
in the framework of general relativity and in a bi-metric theory of gravity proposed 
by Carlip [M]. The present paper and calculations given in [63| refute the incomplete 
understanding of gravitational physics of the light-ray deflection experiments advocated 
in the PPN formalism [65l [75] and in [91]. The principle of equivalence does imply 
the principle of causality for gravitational fleld and demands the fundamental speed of 
gravity to be flnite (equal to the speed of light in general relativity) which was conflrmed 
in the Solar system experiment of starlight's deflection by the moving Jupiter |60j . 

5.7. A Graphical Representation of the Aberration of Gravity 

Two methods of measuring the position of a light-ray deflecting body (gravitational 
lens) in the sky are: 

(1) directly, by observing light (radio waves) emitted by the body itself. This does not 
require any measurement of the gravitational light-ray deflection pattern and relies 
upon propagation properties of the electromagnetic fleld alone; 

(2) indirectly, by measuring the pattern of the gravitational bending of light from stars 
(radio sources) located in the sky around the lens, and extrapolating vectors of 
their gravitational displacement to the common point of intersection which is the 
center of mass of the lens (see Fig. [5]). 

The second method actually traces the lines of the gravitational force exerted on photons 
by the body. The force lines become apparently measurable in the plane of the sky due 
to the gravitational deflection of light. This allows us to map the center-of-gravity of the 
body, regardless of whether the body itself emits light or not. This method of "weak- 
gravitational" lensing is used to localize mass-concentration of dark matter in clusters of 
galaxies [711 El] • In the solar system the weak- lensing method can be effectively used 
to determine the aberration of the gravitational force caused by motion of a light-ray 
deflecting body through measuring the retardation of gravity effect that deflnes position 
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of the body on its orbit and the magnitude of the gravitational time delay of light, 
and comparing it against the retarded position of the body obtained from its direct 
radio observations given in the JPL ephemeris [85] . 

The apparent (optical) and gravitational positions of the light-ray deflecting body 
from which it deflects light by its gravity force, should coincide in the frame where the 
body is static. However, in a moving frame, the apparent position of the body is shifted 
in the plane of the sky, because of the aberration of light, by the angle jS = v/c. The 
gravitational position of the body should be also shifted because of the aberration of 
the gravity force lines by the angle = ev / c = v / Cg, as explained in Fig. [61 If general 
relativity is correct, the aberration of gravity force and that of light must be equal and 
both apparent and gravitational positions of the body should coincide in any frame 
irrespectively of its motion. If the Lorentz invariance of gravity is broken, the apparent 
and gravitational positions of the light-ray deflecting body will coincide in a static frame 
but differentiate in a moving one by the angle 5(3 = (3^ — [3. 

Our interpretation of the aberration of gravity arises naturally because the distance 
r = X — z{s) from observer to the gravitating body can be considered as a spatial part 
of the null vector r" = (r°, r) with r° = Cg(t — s), where t is time of observation, and s is 
retarded time defined by equation (!68l) . The null vector transforms in accordance with 
the group of the Lorentz transformation of gravitational field described by the matrix 
A'^(e) and its direction defines the null direction of the propagation of gravitational 
force from the light-ray deflecting body to the light particle. This null direction of the 
gravity force propagation changes (aberrates) when one goes from one frame to another. 
The magnitude of the aberration of gravity effect is linear with respect to v/c^ like the 
aberration of light is linear with respect to the ratio of v/c. 

This reasoning makes it evident that the measurement of the aberration of gravity 
can be obtained by measuring the deflection of light (radio waves) in a frame that is not 
stationary with respect to the light-ray deflecting body. Our 2002 Jupiter deflection 
experiment [60] was specifically designed to measure the speed of gravity force by 
observing the magnitude of its aberration via excess to the gravitational Shapiro time 
delay caused by retardation in propagation of gravity from Jupiter to radio photons 
from quasar (see text below). Future optical deflection observations by space missions 
SIM [TCy or Gaia [71] of stars near the limb of Jupiter or Saturn, should bring about 
more accurate results. Combining VLBA and Square Kilometer Array [72J to the 
ineterferometer with inter-continental baseline can be also used to conduct this type 
of the gravitational experiments with unparalleled degree of precision |73j . 

The Sun is another moving body with respect to the earth as observed from the 
geocentric frame, but only radio observations of radio sources near the Sun that can 
reach accuracy being sufficient to measure the aberration of gravity effect, are currently 
feasible (see [67] for future prospects). Previous solar bending observations have been 
always analyzed in the barycentric frame (with respect to which the Sun is almost static 
so that the gravity aberrational term is negligibly small) because the main goal was the 
precise measurement of the PPN parameter 7 |102[ 1103] . However, the analysis of the 
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experiments directly in the geocentric frame, where the Sun moves with velocity 30 
km/s, can determine not only the PPN parameter 7, but the gravity aberration term 
proportional to 6a{e), derived below and shown in equation ( IHOl) . This is because the 
Lorentz transformation from the barycentric to geocentric frame probes transformation 
properties of the gravitational force which are controlled by the speed of gravity Cg. 



6. Experimental VLBI Measurement of the Aberration of Gravity 
6.1. The Magnitude of the Aberration 

The formulae for the aberration term can be determined from equations fl68l) and (169|) . 
Let us assume that the impact distance d of the light ray to the massive body is small 
{d ^ r) so that we can expand equation fl69l) in a Taylor series with respect to the 
parameter d/r. Then, we obtain [8| IST] 

r-,.r.^^0{^). (70, 

where d{e) = |^(e)| is the impact distance of the light ray with respect to the retarded 
position z{s) of the body, ^(e) = k x [r{s) x k], r{s) = x — z{s), r{s) = \r{s)\, x is the 
position of observer at the time of observation t, the retarded time s = s{e) is calculated 
from equation fl68l) that, after accounting for definition Cg = c/e, acquires the following 
form 

s = t-e'-^. (71) 

c 

Using equation (jTOl) and neglecting the explicit velocity-dependent terms, we can 
simplify equation for the total angle of the light deflection and reduce it to the 
form 

«(.) = 2(1+,)— Iji. (72) 

where the gravitational deflection of light vector a. = (a*), and M is mass of the body. 

If the speed of gravity Cg = c/e is different from the speed of light c (e 7^ 1) the 
Lorentz-invariance of gravity is broken and the aberration of gravity differs from that of 
starlight. This would produce a small tangential (in the plane of the sky) displacement 

5a.{e) = a(e) — cx (73) 

of star's position, deflected by the body as measured in the moving frame, compared 
with that ex. = ct{e = 1) expected from general relativity. The approximate position 
and motion of the body, however, must be known in order to estimate an accurate value 
for a which is taken as a reference in measuring e = c/cg. For solar system objects, 
exceedingly accurate ephemerides are available [85ii86j, so that a can be predicted with 
an accuracy being sufficient to set a stringent upper limit on the Lorentz violation of 
gravity field and its fundamental speed Cg, which is not limited by the uncertainty in 
position and/or velocity of the body. 
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The aberration of gravity effect is estimated after expansion of a.{e) in a Taylor 
series around e = 1 and making use of equations flTTj) and fl73|) . One has (see Fig. [7]) 

da(e) 



a(e) = a. 



de 



i^-l) + [{e - If] , (74) 

J e=l 

that after substitution to equation ( |73l) and taking the derivative yields 

<5a(e) = ^5/3(e), (75) 

a =2(1 + 7)^, (76) 

e =-, (77) 
r 

5/3(e) =A-/3=(e-l)^, (78) 

c 

where i>t = x (v x fc) is the transversal velocity of the light-ray deflecting body in 
the plane of the sky and all values of the functions entering equations (|751) - (rr5l) are 
taken at the retarded instant of time for e = 1. Clearly, the equation flTHl) compares 
the aberration of gravity, /3e = ^'r/cg, with respect to the aberration of light, f3 = vt/c 
which are both linear with respect to body's velocity v but have different physical origin 
associated with the transformation properties of the gravitational and electromagentic 
field. Measuring presumable tangential displacements Sol in the gravitational deflection 
of starlight sets a stringent limit on the aberration of gravity (3^ and, hence, on a possible 
discrepancy between the speed of gravity Cg and that of light c. 

Equation (175|) was used in our papers [59l [60] to estimate the aberration of the 
gravity force in the light bending experiment with Jupiter as the gravitational lens. We 
analyzed the data in the barycentric frame of the solar system where Jupiter has an 
orbital speed of 4.5 x 10~^ c. The maximal magnitude of the aberration of gravity force 
caused by the motion of Jupiter reaches 5a ~ 10 mas when Jupiter is at the distance 6 
astronomical units (6 AU) from the earth and light passes near Jupiter's limb. For the 
2002 experiment, the quasar passed 3.7' from Jupiter and the deflection of light caused 
by the aberration of gravity force was 0.05 mas. It was measured to an accuracy of 
20% [60j, thus, disproving experimentally the erroneous statements by T. van Flandern 
[95| [96] about the nature of gravity and its speed of propagation. 

Most of the previous radio gravitational bending experiments occurred in early 
October when the Sun passed in the plane of the sky closely in front of the strong quasar 
3C279. The aberration of gravity term as viewed in the Shapiro time delay computed 
in the earth-centered frame is relatively large since the Sun moves with a f © ~ 30 km/s 
in this frame. The estimate of the tangential component in the gravitational bending of 
light caused by the aberration of gravity relative to the aberration of light effect, yields 



1+7^ f A \ [Rq 



2 



=«0 W , (79) 
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where a© = 1.75" is the hght deflection on the solar hmb, A = 1 AU is one astronomical 
unit, Rq = 7 X 10^° cm is radius of the Sun. Substituting in the numbers one obtains 

6a{e) = 37.5(e - 1) mas . (80) 

If we assume the minimal value of the impact parameter that is permitted in radio 
observations by the solar corona as ~ ARq, the maximal value of the aberration of 
gravity effect (see Figure [7]) is 6a ~ 2344 /xarcsec for e = 0. Assuming that the precision 
of VLBI measurement of the quasar's position is 20 /iarcseconds we shall be able to 
measure the aberration of gravity effect in the solar light-ray deflection observations 
with the accuracy approaching, in principle, to 1% as contrasted to 20% in the case of 
the Jovian deflection experiment. 



6.2. Lorentz Transform Between Frames 

The fundamental measurable of VLBI is the difference in the arrival time of radio waves 
from an external source impinging on each telescope (station) in the array. This time 
difference is inferred by the phase difference of the two radio waves. The source of radio 
waves can be a quasar at cosmological distance or a spacecraft in the solar system. The 
former produces plane waves and is stationary in the sky, whereas the latter is in the 
solar system which adds minor complications to the analysis that are not considered 
here. 

Let us consider two frames, the barycentric frame of the solar system and the 
geocentric frame, on which VLBI array analyses are done. As a quasar signal passes by 
a radio- wave deflecting body (Sun, planet) on its way to the array, the major components 
of the observed delay (difference in the arrival time between two stations, as measured 
the phase difference of the two signals) in the barycentric frame are given by equation 

m 

T2 - Ti = ~-K B + H{T, X) + E{T, X, v) , (81) 

R, =X,-Z(5,), (z = l,2), (83) 

Si =T,-^, {^ = l,2). (84) 

where we use capital letters T and X to denote the time and spatial coordinates in 
this frame, Tj {i=l,2) is the time of arrival of radio wave to i-th VLBI stations, K is 
the unit vector toward the radio source, H is the gravitational time delay, E represents 
all additional sources of delay error and is (radio) frequency u dependent (egs. solar 
coronal refraction), Xi = X{Ti) is coordinate of i-th VLBI station at the time Tj, 
with B = X2 — Xi as the basehne between the two VLBI stations, Z{Si) is a retarded 
coordinate of the light-ray deflecting body at the retarded time Si, V{S) = dZ{S)/dS is 
the barycentric velocity of the body taken at the retarded time, and Ri is the distance 
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between the station i and the retarded position of the hght-ray deflecting body. In 
§VIB.5 we will outline the methods that are used to isolate the gravitational delay, H, 
by removing the delay errors E. The solar light-deflection experiments neglect the slow 
motion of the Sun with respect to the barycenter, so that its coordinate Z is considered 
as flxed during the time of the experiment. The static- fleld relativistic correction H to 
the time delay was derived by [50]. Note that the static- fleld approximation can not 
be applied in the deflection experiment of Jupiter in the barycenter frame for it moves 
with respect to this frame rather rapidly; hence, the aberration of gravity effect from 
its velocity in this frame was taken into account in the analysis of the September 2002 
experiment with Jupiter [591 160] . 

The Lorentz transformation of the electromagnetic and gravitational flelds to a 
frame moving with velocity v with respect to the barycentric frame, using equation 
( 167|) . gives the following components for the gravitational delay in the moving frame 

^2 - ti = --k b + hit, x) + e(t, X, u') , (85) 

n =x,-z{si), (^ = 1,2), (87) 
Si =ti--, {t = l,2). (88) 

where t and x denote the time and spatial coordinates in this moving frame, (i=l,2) 
is the time of arrival of radio wave to i-th VLBI stations, k is the unit vector toward 
the radio source, h is the gravitational time delay, e represents all additional sources of 
delay error and is frequency z/' dependent (notice that v' ^ v because of the Lorentz 
transformation), Xi = x{ti) is coordinate of i-th VLBI station at the time tj, with 
b = X2 — xi as the baseline between the two VLBI stations. In converting to the 
moving frame, the speed of gravity, Cg, appears as a consequence of transformation of 
the gravitational fleld variables in two places: flrst, in the aberrational term k ■ v/cg in 
front of the logarithm; and second, in the retarded term r^, the distance between the 
station i and the light-ray deflecting body taken at the retarded time Sj. 

The general form of the time delay equation which takes into account not only 
linear but quadratic and high-order velocity-dependent corrections as well, has been 
derived in |8l [51] and [83] under assumption that Cg = c. 

6. 3. Measuring the PPN Parameter 7 

The gravitational time delay is proportional to (1 + 7), where the parameter 7 is a 
static measure of space curvature and is unity in general relativity and zero in the 
Newtonian gravity [64J. The quantity (7 — 1) measures the degree to which gravity is 
not a pure geometric phenomenon and other gravity-generating scalar flelds are involved 
[68] . These flelds might dominate in the early Universe, but would have now weakened 
so as to produce tiny, but presumably detectable effects [771 168]. 
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Although all VLBI experiments are done on or near the earth, most astrometric 
quantities, including 7, are more conveniently analyzed in the solar system barycentric 
reference frame (T, X) because most velocity- dependent relativistic corrections are 
zero in this frame. Hence, VLBI reduction systems use this frame [88l [89l [901 l9T] 
where the gravitational time delay is given by equation (l82ll with V = 0, since this 
velocity is very small. Thus, in the solar light-ray deflection experiments, we should 
analyze the measurements of the gravitational delay as a function of the Sun-quasar 
angular separation in the barycentric frame to determine 7. The estimated accuracy in 
measuring 7 — 1 each observing day for this kind of experiments with 3C279 that passes 
near the Sun each October, is given in Table 1 and is expected to be about 7x 10~^. The 
best previous measurement of 7 = 2.3 x 10~^ was achieved during the conjunction of the 
Sun and Cassini spacecraft in 2002 using the gravitational-induced frequency change of 
the spacecraft transmitters as Cassini went behind the Sun on its way to Saturn [97] . 

6.4- Measuring the Aberration of Gravity and Its Fundamental Speed 

In the 2002 Jupiter deflection experiment, the aberration of gravity in the retarded 
position of Jupiter was measured in the barycentric frame by making use of equations 
(!8T1) - (|8^ because Jupiter moves with respect to the barycenter of the solar system [60] . 
But, the barycentric velocity of the Sun is rather small, so its gravity- field aberration 
in the barycentric frame is not easy detectable. However, the aberration of the solar 
gravity field, is measurable from the earth-center frame since the deflecting object, the 
Sun, is moving about 30 km/s with respect to this frame and its gravitational field 
contains the gravimagnetic component that drags photons in the direction of its motion 

The experimental observable is the measured delay difference of the electromagnetic 
wave between two stations. This delay (technically, it is a phase difference, but can be 
converted into delay) is a scalar function of time and baseline. In the absence of gravity 
field this delay is Lorentz invariant, that is [631 1211 ES] 



T2-Ti + -K-B = t2-ti + -k-b. (89) 
c c 

Therefore, the gravitational part of the time delay difference must be also independent 
of the observing frame used in the analysis (because it is a scalar), 

H{T,X) = h{t,x) , (90) 

where the conversion of all relevant quantities from the barycentric to the geocentric 
frame must be done precisely. 

In the following description of the conversion of the reference frames, we will denote 
the barycentric coordinates of the Sun Zq, those of the geocenter Z^, and the geocentric 
coordinates and velocity of the Sun are Zq = Zq — Z^ and Vq = Vq — Vq respectively. 
Relativistic corrections to these definitions are quadratic with respect to v/c [6811891190] 
and can be neglected in current practical application of equation (1901) . The barycentric 
coordinates and velocity of the Sun and the geocenter are also known as functions of the 
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barycentric time T from accurate solar system ephemerides [55] . The barycentric 
time T is related to the geocentric time t measured on the earth by a time transformation 
in which the difference is quadratic with respect to v/c terms |100l llOlj . and, thus, 
negligibly small in dealing with equation (!90|) . The Lorentz transformations needed for 
the frame conversions in equation ( |90ll cirG clS follows. 

First, the two unit vectors K and k, characterizing the direction of the propagation 
of light in two different frames, are connected by the Lorentz transformation equation 
( l39l) for electromagnetic field. In the linear approximation this transformation is reduced 
to the classic expression for the aberration of starlight 

k = K--Kx{V^xK) + (^] . (91) 



c ■ " ■ V^^/ 
Barycentric earth's velocity and the astrometric coordinates of the vector K for 
3C279 and the other calibrators are accurately known, so the coordinates of the unit 
vector k in the geocentric frame can be easily calculated from equation (1911) . and vice 
versa. 

Second, the retarded geocentric coordinate of the moving Sun, ZQ{si), entering 
equation (186|) . can be represented in the first approximation as a difference between the 
barycentric coordinates of the Sun, Zq, and the geocenter, Z^, 

= Zeis,) - Zeis,) + O (^] , (92) 



where the retarded time 



Si = ti - —\xiti) - ZQisi)\ . (93) 

Next, the coordinates of the VLBI station in the moving and static frames are related 
by equation 

x, = X,(t,)-Ze(t.) + 0(^^^ , (94) 

taken at the time of observation tj. The retardation in the barycentric coordinate of the 
center of mass of the Sun Zeis,) does not affect its value significantly because it moves 
very slowly around the barycenter of the solar system; hence equalities Zeis,) = Zeiti), 
Veisi) = Veiti) can be used, although the precise correction for the retardation is made 
in the analysis software to eliminate any possible error. On the other hand, the orbital 
motion of the earth around the barycenter of the solar system is significant and the effect 
of the retardation of the force of gravity (as it propagates in the geocentric frame from 
the Sun to the earth) in the coordinates of the geocenter must be taken into account. 
It yields for the retarded barycentric coordinate of the geocenter, 

Zeis,) = ZM + Veit,)is, -ti)+0 [isi - t,)'] (95) 



1 

= Ze(t.) - -Veiti)\xi - ZQisi)\ +0(- 

Cg 

where the retardation of gravity correction proportional to Si — ti, is linear with respect 
to the orbital velocity of the earth and has the same order of magnitude as that for 
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the aberration of light given in equation (MT\i . if Cg = c. Finally, transformation of the 
spatial coordinates from equation fl87|) to fl83l) becomes 



Introducing unit vectors = r^/rj and Ni = Ri/Ri we can re- write equation (19^ to 
the following form 



Note, that the transformation of the retarded coordinates of the earth and the the 
Sun relates exclusively to the transformation of the gravitational field variables, while 
the transformation of electromagnetic field (radio wave from the quasar) is associated 
with the change fPTl) in the direction of its propagation. Thus, equation (jHT]) describes 
the aberration of gravity, while equation (19T|) is the aberration of starlight (see Fig. 
[6]). These aberrations could be different if the fundamental speed of gravity Cg were not 
equal to the speed of light c. This difference would lead to the gravitational deflection of 
light coming from a quasar, different than the prediction of general relativity as shown 
in Fig. [71 

The velocity- dependent term in front of the logarithmic function in equation (!82l) 
is proportional to the radial velocity of the Sun with respect to the barycenter, K ■ Vq. 
This velocity is too small and its effect on the time delay H in the barycentric frame 
is not currently detectable. On the other hand, the radial velocity-dependent term in 
front of the logarithm in equation (!86|) is k ■ {Vq — V^). It includes the orbital velocity 
of the geocenter and, in principle, could be measured. However, the orbital velocity 
of the geocenter is projected on the direction to the source of light k, and reduces its 
effect on the time delay to the size of about 1 yuarcsecond which is somewhat too small 
for current VLBI technology. Therefore, we shall neglect the velocity-dependent terms 
in front of the logarithmic functions. 

After the Lorentz transforms are made from the barycenter to the earth-center 
frame, we can measure the fundamental speed of gravity Cg by solving equations 
fl85|) - fl88l) for the gravity Lorentz-invariance violating parameter e. An estimate of its 
measured accuracy from the solar deflection experiment is a few percent and depends 
on how accurately the coronal refraction is removed. 

We emphasize that the above transformations must be made as accurately as 
possible, given the desired accuracy of the observations. For example, the unperturbed 
light particle moves in vacuum with constant velocity c in any frame. Since time can 
be measured presently with much better accuracy than space intervals, this led to 
abandoning the measurement of the speed of light c in vacuum more and more precisely. 
Thus, the numerical value of the speed of light has been fixed c = 299792458 m-s~^, and 
it is this value which is used to measure distances in the solar system. The most accurate 
barycentric coordinates of the solar system bodies - Sun, planets, and their satellites 
- are normalized to this numerical value of the speed of light c and are tabulated in 




(96) 




(97) 
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the solar system ephemeris [HSj- We should adopt this metric system of units and 
consider the speed of light c as known parameter that is not subject to measurement. 
This point is of paramount importance for unambiguous physical interpretations of 
gravitational experiments in the solar system and beyond [62]. 
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Figure 1. Gravitational field of a moving body is a retarded solution of the gravity- 
field wave equation (fT3|) . (fT4|) . In general relativity, gravitational field originates at 
the body location at each instant of time s on its world line at the point z(s), and 
progresses on the hypersurface of a null cone, attached by its vertex to the point z{s), 
from the past to the future. Directions of the propagation restricting the domain of 
the causal influence of the gravitational field of the body located at the point z{s), are 
shown by arrows. The conventional interpretation of this drawing is that the moving 
body "emits" the gravitational field at time s, and observer measures the gravitational 
field at time t, when the body is located at the retarded position z{s) on its orbit at 
the retarded time s = t — r/cg, which is a null characteristic of the gravitational field 
as follows from equation (jlSp . For simplicity, the picture shows only one null cone 
of gravity field. Effectively, the entire sequence of null cones must be shown for each 
instant of time corresponding different positions of the body on its world line. 
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Observer 
51V, 




Light-ray-deflecting body's worldline 



Figure 2. A photon is emitted by a source of liglit (star/quasar) and propagates along 
a hypersurface of light null cone (solid red line) toward observer. Gravitational field of 
a light-ray deflecting body propagates along the gravity null cone (dashed green lines) 
and interacts with the photon at the points of intersection of the light and gravity null 
cones (points 1 through 5). Gravitational light-ray deflection angle and time delay are 
integrals along the light null cone but their numerical values are determined by the 
retarded position of the body on the gravity null cone connecting the body and the 
photon as follows from equations (p7)) and Observer measures the time delay and 
the deflection angle of the photon at time t corresponding to point 5. At this time 
the time delay and the deflection angle are determined by the body position taken on 
the "last" gravity null cone at point E. Notice that the point of the closest approach 
of light ray and/or the time of the clossts approach of the photon to the body are 
irrelevant to the process of interaction of light with gravitational field of a moving 
body, as follows from the result of integration of the light-ray propagation equation. 
This picture illustrates the case when the gravity and light null cones are formally 
different (cg 7^ c). In general relativity the gravity null cone coincides with the light 
null cone and the photon, after it passes the body, moves in the gravitational field of 
the body which is practically " frozen" for this photon ^51, because the speed of gravity 
is the same as the speed of light. 
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Figure 3. Gravitational field of a moving body affects only the particles lying on 
the hypersurface of the future gravity null cone due to the gravity causal nature. A 
photon emited by a source of light (star, quasar) at time to arrives to observer at time 
t along a null direction of the past light cone with the observer located at its vertex. 
As explained in section 5 and in caption to Fig. [21 the photon detected at time t, is 
deflected by planet's gravity force from the planet's retarded position taken at time 
s = t — r/cg that is a null characteristic of the retarded solution of the gravity- field 
wave equations (fT^ . This effect of the retardation of gravity can be observed 

by measuring the amount and direction of the gravitational deflection of starlight by 
a moving body. In general relativity Cg = c, hence, the future gravity null cone of the 
body and the past light cone of the observer must coincide along the null direction 
which is simultaneously characteristic of the Maxwell and Einstein equations. This bi- 
characteristic nature of general relativity led to the confusing statements in literature 
(see [65l [75] and references therein) that the Jovian light-ray deflection experiment 
measured the speed of light coming from the quasar. In fact, the light was used to 
measure the amount and direction of its gravitational deflection from which the orbital 
position of the Jupiter was derived as predicted by the equations (|67|) . (|68|) [60j . 
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(A) (B) 

Figure 4. Light-ray deflection by a static (A) and moving (B) body. In case (A) tlie 
spatial distance between the body and observer does not change as light propagates. 
Thus, measuring the deflection of light does not allow us to determine experimentally 
whether the gravity force of the body acts with retardation from position (1), or 
instantaneously from position (2). In case (B) the spatial distance between the body 
and the observer varies as the photon travels toward observer. The retarded interaction 
of gravity with light becomes apparent since measuring the angle of the gravitational 
deflection of light allows us to distinguish between positions (1) and (2) of the body 
on its world-line making the causal structure of the gravity cone clearly discernible 
and measurable. This kind of gravitational experiments works perfect even if the 
body moves uniformly with constant velocity with respect to observer so that the 
gravitational radiation is not emited. 
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Figure 5. Deflection of starlight caused by spherical gravitational field of a body. 
Arrows show the angular displacement of stars from their undeflected positions and 
map the lines of the gravity force deflecting the starlight. The point of intersection 
of these lines deflnes the body's center of gravity in the sky. Apparent position of 
the body visible in optics/radio and that found from the gravitational deflection of 
starlight coincide in the static frame. 
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Gravitational position of the lens in a moving frame 




Position of the lens in the static frame 
The aberration of gravity angle p 



fl^ ' Optical position of the lens in a moving frame 



The aberration of light angle (3 



Figure 6. Deflection of light in a moving frame. The apparent optical positions of 
the stars and the light-ray deflecting body are shifted by the aberration of light from 
their static-frame positions. If the Lorentz invariance of gravity is broken, the body's 
center of gravity, determined from the gravitational deflection of light, will not coincide 
with its optical position. Light-ray deflection experiment conducted in a moving frame 
measures the difference 6(3 = Pe — defined in equation ([78|l . In general relativity 
(5/3 = 0. 
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Deflected Position of 3C279 
in general relativity 



Gravitational position of the sun 
in the geocentric frame 



5a 

Deflected Position of 3C279 in 
' 7 the theory of gravity with broken Lorentz invariance 

a / 

Undeflected Position of 3C279 



Optical position of the sun 
in the geocentric frame 



Figure 7. The difference between the aberration angles for gravity force and light in 
the plane of the sky in the geocentric frame. The angular separation of 3C279 from the 
Sun is 6 — (i/(lAU). The apparent position of the quasar is shifted in the sky by the 
force of solar gravity at the angle a if Cg = c. In case oi Cg ^ c the quasar is displaced 
by gravimagnetic force in the direction of orbital motion of the earth at a small angle 
5a given by equation ((79|) . It corresponds to the angle 5(3 = [e ~ 1){vq/c) between 
the directions to the gravitational and optical positions of the Sun. The relationship 
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TABLE 1: ESTIMATED EXPERIMENTAL SENSITIVITY TO 7 and ABERRATION 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Date 


(Sol Rad) 


(mas) 


(mas) 


(mas) 


(mas) 


(10-5) 

V / 


(10-2) 


October 1 


28.9 


35 


0.0 


0.039 


0.015 






October 5 


10.0 


176 


0.1 


0.375 


0.015 


21 


4.0 


October 6 


6.3 


279 


0.3 


0.945 


0.015 


12 


1.5 


October 7 


2.5 


704 


9.1 


6.000 


0.090 


33 


1.5 


October 9 


5.2 


338 


2.1 


1.386 


0.025 


13 


1.8 


October 10 


9.0 


196 


0.2 


0.363 


0.015 


13 


4.1 


October 11 


12.8 


137 


0.1 


0.229 


0.015 


19 


6.6 


October 18 


39.6 


40 


0.0 


0.018 


0.015 






All days 












7 


1.0 



(1) = Mean observing time 

(2) = Sun-3C279 separation 

(3) = Gravitational bending of light 

(4) = Estimate coronal bending at 23 GHz 

(5) = Gravitational aberration of Sim 

(6) = Expect positional sensitivity; 0.015 + 1% of coronal bending 

(7) = Accuracy of 7 — 1 

(8) = Accuracy of aberration of gravity e — 1 



